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A nonlinear stochastic heat equation: Holder continuity 
and smoothness of the density of the solution 

Yaozhong Hu* David Nualart^ and Jian Song 

Abstract 



In this paper, we establish a version of the Feynman-Kac formula for multidimen- 
sional stochastic heat equation driven by a general semimartingale. This Feynman-Kac 
ly-^ I formula is then applied to study some nonlinear stochastic heat equations driven by 

CN ■ nonhomogenous Gaussian noise: First, it is obtained an explicit expression for the 

Malliavin derivatives of the solutions. Based on the representation we obtain the 
I smooth property of the density of the law of the solution. On the other hand, we also 

Ph ■ obtain the Holder continuity of the solutions. 

■ 

1 Introduction 

^ . In this paper we consider the following nonlinear stochastic heat equation: 

— ■ \-^ = -Au + b{u) + a{u)W{t,x), t>0, xeR" 

[u{0,x) = Uo{x) , 

O ■ ^ 2 

where A = J2i=i ^ the Laplace operator, b and a are globally Lipschitz continuous 
functions, and is a zero mean Gaussian random field, which is a Brownian motion in the 
time variable and it has a nonhomogeneous spatial covariance with density q{x, y) (see (12. ip 

^ ■ for the precise definition). Here W{t,x) denotes the generalized random field — — — . 

5^. The case of an homogeneous covariance kernel q{x,y) = q{x — y) has been studied 

in the seminal paper by Dalang [3]. In this case, the existence, uniqueness and Holder 
continuity of u{t,x) with respect to both parameters t and x is obtained in under 
integrability conditions on the spectral measure /i of the noise. We extend these results to 
the nonhomogeneous case in Section 4. 

On the other hand, using the techniques of Malliavin calculus, and assuming suitable 
nondegeneracy conditions, one can show that for a fixed {t,x), t > 0, the random variable 
u{t,x), solution to (11.11) . has an absolutely continuous probability law and the density is 
smooth. The results that have been obtained so far along this direction can be summarized 
as follows. 
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(i) In [13] Pardoux and Zhang considered Equation fll.ip when x is in the interval (0, 1) 
with Dirichlet boundary conditions, assuming that W is a. space-time white noise. 
In this case, if the coefficients are Lipschitz, then u{t,x) has an absolutely continuous 
distribution for any t > provided a{uo{xo)) ^ for some x^ e (0, 1). The smoothness 
of the density in this framework was proved by Mueller and Nualart in |TD], assuming 
that the coefficients are infinitely differentiable with bounded derivatives. On the 
other hand, under the stronger nondegeneracy condition > c > 0, and smooth 
coefficients. Bally and Pardoux [1] proved that the law of any vector of the form 
(^(t, xi), . . . < xi < ■ • • < < 1, t > 0, has an infinitely differentiable 
density, assuming Neumann boundary conditions on (0, 1). 

(ii) For the d-dimensional heat equation with an homogeneous spatial covariance, Nu- 
alart and Quer-Sardanyons have provided sufficient conditions for the existence and 
smoothness of the density of uit^ x) for t > and x G M'^, assuming |cr(x)| > c > 0, in 
the paper |12j (see also [1]). 

An open problem for the stochastic heat equation with colored spatial covariance is to 
derive the existence and smoothness of the density under a nondegeneracy condition of the 
form (t(mo(xo)) 7^ for some xq G M'^. The main purpose of this paper is obtain new results 
in this direction. To prove such results we need to show that the norm of the Malliavin 
derivative of the solution \\Dsu{t, x)\\1i^ds is either strictly positive almost surely (for the 
absolute continuity) or it has negative moments of all orders (for the smoothness of the 
density), where "H is the Hilbert space associated with the spatial covariance. 

We develop a new approach to prove these results based on the Feynman-Kac represen- 
tation for the solution to the heat equation with multiplicative noise driven by a general 
continuous semimartingale. The main idea is to express \\Dsu{t, x)\\'^ as the norm in L^(M'^) 
of a function Vs^^{t,x) given by Vs^^{t,x) = J^ac{(,,y)Ds^yu{t,x)dy, where c is the square 
root of the kernel q as an operator. Then for any fixed (s,^), Vs^^{t,x) satisfies the linear 
stochastic heat equation with random coefficients 

= - AK,c + + a'{u)V,,^W{t, x), t>s,xeR'', (1.2) 

with initial condition Vs,^{s,x) = c{^,x)a{u{s,x)). 

In order to establish a Feynman-Kac representation for the solution to Equation (11. 2p 
we need to assume that the covariance kernel q{x,y) is non-singular and this implies the 

dWi 

existence of a random field Wi{t,x) such that W{t,x) = — — {t,x). Then, Equation (II. 2p 

at 

is a particular case of a more general stochastic heat equation of the form 

dV 1 dF 

^{t,x) = -AVit,x) + V^{t,x), (1.3) 

where {F{t,x),t > 0,x G M'^} is a continuous semimartingale in the sense of Kunita [9], 
with local characteristic b(t,x) = b'{u(t,x)) and a(t,x,y) = a'{u(t,x))a'{u(t,y))q{x,y). In 
Section 3 (see Theorem 13. ip we derive a Feynman-Kac formula for the solution of (II. 3p 
assuming that the functions b and a are bounded by C(l + |a;|^), and C(l + + \y\^) for 
some < /3 < 2. This result has its own interest. The proof is based on a generalized Ito 
formula proved in [9]. 
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There have been other papers on the Feynman-Kac formula for the stochastic heat 
equation. We can mention the recent works [6] and [7] on the stochastic heat equation 
driven by fractional white noise. We refer to the references in these papers for related 
works. 

In Section 4 we show the existence and uniqueness of a solution for the general stochastic 
heat equation (11.11) with a nonhomogeneous spatial covariance and we deduce the Holder 
continuity of the solution. This result is an extension of the results proved in [15]. Finally, 
in Section 5, assuming that the covariance kernel is continuous and under a nondegeneracy 
condition of the form q{xo,Xo) > and cr(Mo(a;o)) 7^ for some Xq G M, we estabhsh the 
absolute continuity of the law of the solution and the smoothness of the density if the 
coefficients are smooth. 

To simplify the presentation we have assumed that the functions b and a depend only on 
the variable u. All the results of this paper could be extended without difficulty to the case 
of coefficients b{t, x, u) and x, u) such that they are Lipschitz and with linear growth in 
M, uniformly in (t, x) G [0, T] x R'^ for any T > 0. In this case, the nondegeneracy condition 
would be (t(0, Xo, Mol^^o)) 7^ 0, for some Xq G W'-. 

The results of this paper can be extended to the stochastic heat equation on an open and 
bounded set A C M*^, with Dirichlet boundary conditions. In this case, the Feynman-Kac 
formula involves a d-dimensional Brownian motion starting form a point x G A, and killed 
when it leaves the set A. On the other hand, the existence and smoothness of the density 
have been deduced, applying techniques of Malliavin calculus, for stochastic differential 
equations of the form Lu = h{u) + a{u)W , where L is a differential operator more general 
than 9f — |A (see, for instance, [8l[T6] where L is a pseudodifferential operator and |T2] where 
L is a general parabolic or hyperbolic operator). In all these examples, one assumes that 
a is bounded away from the origin. Our approach to handle a nondegeneracy of the form 
cr(0, Xo, Mo(xo)) only works if a Feynman-Kac representation is available for the corresponding 
stochastic linear equation satisfied by the Malliavin derivative. This happens, for instance, 
for parabolic operators of the form L = dt — bidx^ — \ Yliij (^i,j^Xi,xj- The methodology 
developed in this paper could be extended to these operators, replacing the Brownian motion 
by the diffusion process with generator bid^^ — \ Ylii j ^i,j9x^,xj- 

2 Preliminaries 
2.1 Malliavin calculus 

Let {Q,J-',P) be a complete probability space. Consider a family of zero mean Gaussian 
random variables W = {Wt{ip),ip G C^{R'^),t > 0}, where C~(R'^) denotes the space of 
infinitely differentiable functions on M"^ with compact support, with covariance 



where g is a nonnegative definite and locally integrable function. 

Let "H be the Hilbert space defined as the completion of C^{M.'^) by the inner product 




(2.1) 
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The mapping l[o,t]V ^ Wt{ip) can be extended to a linear isometry between T-Loo '■ = 
L^([0, oo); "H) and the space spanned by W. Then {W{h),h G Tioo} is an isonormal 
Gaussian process associated with the Hilbert space Tioo- 

We will denote by D the derivative operator in the sense of Malliavin calculus. That is, 
if F is a smooth and cylindrical random variable of the form 

F = f{W{h),...,W{hn)), 

hi G Tioo, f G C^(R") (/ and all its partial derivatives have polynomial growth), then DF 
is the Hoo-valued random variable defined by 

DF = J2§^,iW{h),...,W{K))h,. 

The operator D is closable from L^(fi) into LP'{VL] T-Loo) and we define the Sobolev space D^'^ 
as the closure of the space of smooth and cylindrical random variables under the norm 

ll^^lli,2 = ^Je{f^) + e{\\df\\Ij. 

We denote by 5 the adjoint of the derivative operator, given by duality formula 

E{Siu)F) = E{{DF,u)^J, (2.2) 

for any F G D^'^ and any element u G "Hoo) in the domain of 6. The operator S is also 

called the Skorohod integral. The higher Malliavin derivatives can be defined in similar way 
and we can define ©'^'^ for any integer k > 1 and real number p > 1. Set = P| D*^'^. 

k>l ,p>2 

To obtain the existence and smoothness of the density, we make use of the following criteria. 

Theorem 2.1 Let F : f2 — )■ M 6e a random variable. If F G D^'^ and \\DF\\y^^ > almost 
surely, then the probability law of F is absolutely continuous with respect to the Lebesgue 
measure. Moreover, z/F G and E [||F)F||^^] < oo for all p > 1, then the density of F 
is infinitely differentiable. 

For the proof of this result and a detailed presentation of the Malliavin calculus we refer 
to [11] and the references therein. 

2.2 Generalized Ito formula 

In this section we introduce some preliminaries on continuous semimartingales depending 
on a parameter and the corresponding generalized Ito formula. We refer to [9] for more 
details. 

Fix a time interval [0, T], a complete probability space (fi, J-", P) and a filtration {J-'j, < 
t < T} satisfying the usual conditions (increasing, right-continuous, and J-q contains all the 
null sets). Let {F(t,x),0 < t < T, x G 0} be a family of real valued processes with 
parameter a; G O, where O is a domain in W^. We can regard it as random field with 
double parameters x and t. If F{t,x) is m-times continuously differentiable with respect to 
X a.s. for any t, it can be regarded as stochastic process with values in or a C™-process. 
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Here we denote by C"^ = C™(0,M) the set of all real valued functions on O which are m 
times continuously differentiable. If furthermore, for each multi- index a G {1, . . . , c?}^ with 
\a\ = k < m, {D"F{t,x),x G 0} is a family of continuous semimartingales, then F(t,x) is 
called a C^-semimartingale. Here we have used the notation = 75 — — . 

We denote by C^'^ the set of all functions on a : [0, T] x O x O — M such that the 
partial derivatives ^(t, x, y), ■§f-i't, x, y) and gf.g^. {t, x, y) exist for any 1 < i, j < d and are 
continuous in {x,y), and for any compact set C O and I < i, j < d 

We also denote the set of all functions on 6 : [0, T] x O — t- M which are continuously 
differentiable in x, and for any compact set C O and any 1 < i < d 

^ f db \ 

sup \h{t,x) \ + |- — ] dt < 00. 



x&K 



dxi 



Let {F{t, x), X G 0} be a family of continuous semimartingales decomposed as x) = 
M(t, x) +-B(t, x), where M(t, x) is a continuous local martingale and B{t, x) is a continuous 
process of bounded variation. Let A{t,x,y) be the joint quadratic variation of M(t,x) 
and M{t,y) and assume that A{t,x,y) = a{s, x,y)ds and B{t,x) = Jl^b(s,x)ds, where 
a(t,x,y) and b(t,x) are predictable processes. Then {a(t,x,y),b(t,x)) is called the local 
characteristic of the family of semimartingales {F{t,x),x G O}. Following Section 3.2 of 
[9], we say that the local characteristic (a, b) belongs to the class 5^'° if a(t, x, y) and b(t, x) 
are predictable processes with values in C^'^ and C^, respectively. 

Now let {F(t,x),x G 0} be a continuous semimartingale with local characteristic {a,b). 
Let {ft, < t < T} be a predictable process with values in O satisfying 

/ a{s, fs, fs)ds < 00, / \b{s, fs)\ds < 00 a.s. (2.3) 
Jo Jo 

Then, the Ito stochastic integral of ft based on the kernel F{dt, x) is defined as the following 
limit in probability if it exists 

/ Fids, /,) = lim V{F(t,+i A t, ft,M) - F{h A t, /^^.^J}, 
Jo lAKo^ 

where A = {0 = to < ■ " " < = T}, and |A| = maxi<j<„(tj — ti_i). 

The joint quadratic variation of the Ito integrals F{ds, fs) and F{ds, Qs) satisfies 

F{ds,fs), j F{ds,gs))t= a{s, f^, gs)ds. 
Jo Jo 

The following is the generalized Ito formula (see Theorem 3.3.1 in ^). 

Theorem 2.2 (Generalized Ito formula) Let {F[t,x),x G 0} be a continuous C'^ -process 
and a continuous -semimartingale with local characteristic belonging to the class B^'^ and 
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let {Xt,0 <t< T} be a continuous semimartingale with values in O. Then {F(t,Xt),0 < 
t <T} is a continuous semimartingale and satisfies 

* dF 



F{t,Xt)=F{0,Xo)+ / F{dr, Xr) + y2 / j^{r,Xr)dX, 

Jo j^]^ ^0 (y^i 



for any t G [0, T]. 



2,7 — 1 i — 1 



3 Feynman-Kac formula 

In this section we establish a general Feynman-Kac formula for the dimensional heat 
equation driven by a continuous semimartingale. Suppose that F = {F{t, x),0 < t < T,x E 
M.'^} is a continuous semimartingale with local characteristic {a,b). We are going to impose 
the following condition. 

(HI) Assume that a{t,x,y) and b{t,x) are continuous and satisfy 

\a{t,x,y)\<C{l + \xf+\yf), (3.1) 
\b{t,x)\<Cil + \xf), (3.2) 

for t G [0,r], with 0<(3<2. 

Consider the stochastic heat equation 

dV , , 1, OF 

r(x,0) = h{x). 



— {t, x) = -AVit, x) + V—{t, x) ^3 3^ 



An adapted random field {V{t, x),0 < t < T,x E M'^} is called a mild solution to the above 
equation if V(t, x) satisfies the following integral equation 

V{t,x)= [ pt{x - z)h{z)dz + [ ( [ pt-r{x ~ z)V{r,z)F{dr,z)] dz, (3.4) 



where pt(x) = (27rt) 2 exp( — |xp/2t). 

Theorem 3.1 (Feynman-Kac Formula) Let h{x) be continuous and with polynomial growth. 
Then the process 

V{t, x) = (^{x + Bt) exp (^j^ F{dr, x + Bt - 5,) - ^ ^ a{r, x + Bt - Br)dr^ ^ , 

(3.5) 

where B is a d-dimensional standard Brownian motion independent of F , E^ denotes the 
mathematical expectation with respect to B, and d{t,x) = a{t,x,x), is a mild solution to 
Equation IfST^) . 
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Proof. 

We divide the proof into three steps. 

Step 1. First we show that the process (13.51) is well defined. In the sequel we denote by E 
the mathematical expectation in the probability space where F is defined, and denotes 
the expectation with respect to the independent Brownian motion B. Set 

Yt= [ F{dr, x + Bt-Br)- \ [ a{r, x + Bt - B,)dr. 
Jo ^ Jo 

Notice that j^F{dr,x + Bt — Br) is a well defined Ito stochastic integral, because the 
process {Bt — B^., < r < t} is independent of the semimartingale F, and conditions (12. 3p 
are satisfied. Then V{t,x) = E^ {h{x + Bt) exp(y()). We claim that this expectation exists 
and V{t,x) satisfies the following condition for any x G M'^ and p >1, 

sup E\V{t,x)\P < Kiexp (K2\xf) , (3.6) 

0<t<T 

where the constants Ki and K2 depend on p and T. In particular, this implies that the 
stochastic integral in (13.41) is well defined. We can write 

E\V{t,x)\'P < {E^\h{x + Btl'^^ EE^exp{2pYt))K 
Let us denote by M{t,x) the martingale part of F{t,x). Then we make the decomposition 



Y, = Yt^'' + Y^\ 



where 



and 



r/i) = / M{dr, x + Bt- Br) - p I a(r, x + Bt - Br)dr. 



b{r, x + Bt~Br)+ [p-^] a(r, x + Bt ~ Br 



dr. 



Using conditions (13. ip and (13. 2 p and taking into account that /3 < 2, we obtain for all 

te[o,T] 

E^ exp{2pY}^^) < E^ (^exp (^C {1 + \x + Bt - Br\'^)dr^^ 
< Kiexp{K2\xf). 

On the other hand, taking into account that exp{2pYt^^) is a martingale, we can write 

EE^ exp{2pY}^^) = E^ E exp{2pY}^^) < 1, 
which completes the proof of (13.60 . 

Step 2. We now show that the process (13. 5p is a solution to Equation (13. 3p under some 
additional regularity assumptions on the semimartingale F{t,x). Suppose that F(t,x) is 
a C^-semimartingale, such that the local characteristic (a, 6) satisfies (13.11) and (13.21) . We 
also assume that the functions D"Z)°a(t, x, y) satisfy the estimate (13.10 for all multi-index 



a with 1 < |a| < 2, and the functions D'^b{t,x) and D"d{t,x) satisfy the estimate fl3.2p for 
all multi- index a with 1 < |a| < 2. Clearly this implies that the local characteristic belongs 
to the class 5^'°. Suppose also that the functions D^h(t,x) have polynomial growth for all 
multi- index a with 1 < |a| < 2. 
For fixed x, let 

rt 1 /■*_ 

?/)= / F{dr,x + y — Bj.) — - / a{r,x + y — Br)dr. 
Jo 2 Jq 

According to Theorem 3.3.3 in |9j, ^{t,y) is a C^-semimartingale with local characteristic 
belonging to the class B^'^. We can apply the generalized Ito formula (12. 4p to the process 
Yt = $(t, Bt), and we obtain 



1 r /"* OF 

dYt = F(dt,x) - -d(t,x)dt + y^ / —{dr,x + Bt-Br 

2 ^ Uo dxi 



dBl 



1 

2 ^ 

i=l 

1 



■t=i 



* 9a 
Oxi 



r, X + — Br)dr 



1 r /"* 



'^•^i 



dt. 



The terms (Jg ^{dr, x), B'^)t = (/q ^{dr, x), B'^)t vanish since M and -B are independent. 
The quadratic variation of the semimartingale Y is given by 



d{Y)t = d(t, x)dt + Y \ T—idr, x + Bt - Br 
^-^ axi 



OF 



1 /■* dd 



t^-^i 



2 Jo ^a; 



— (r,x + fit-fi,)f^H (3-7) 



Consider the process Z{t,x) = h{x + Bt)e^^ . Applying Ito's formula to h[x + Bt)e^^ 
yields 



Z{t,x) =h{x)+ [ Z{s,x)dYs + y] [ 
Jo Jo 

1 



9/i 



X + Bs)e^"dB 



dx] 



1=1 



'^•^i 



1 /■* ^ 

X + 5,)e^^cis + - / \/(s,x)d(r), 

2 ./n 



(3.^ 



We claim that the stochastic integrals with respect to B^ in the above expression have zero 
expectation with respect to B. This is a consequence of the following properties 



E''Z{t,xfE'^ 



E 



i=l 

d .t 



dF 





{dr, x + Bt- Br) 



/ E'^Zit.xfE'' y / ^(r,x + Bt- Br)dr 
Jo j^Jo dxi 



dt < 00, 



dt < 00, 



(3.9) 



(3.10) 
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and 



e^^^'ds < oo. 



(3.11^ 



These properties follow from our additional assumptions. For instance, to show (13. 9p for 
the martingale component of F, we take the expectation in the probability space where F 
is defined and we use the fact that for any p > 2 



E 



< CpE 



t o2 



y , 



[r, x + Bt- Br, x + Bt- Br)dr 



<CE (1 + \Bt - Br\^)ds < oo. 
Jo 

Then, taking the expectation with respect to B in (13. 8p yields 



V{t,x) = h{x) + / V{s,x)F{ds,x) 



d 



i=l 



V{s, X 



( n d^F 1 /"'^ d^a 

I a^^^"^' x + Bs-Br)--j^ -^{r,x + Bs- Br)dr 



+ 



+ 



dF 1 da 1^1 

/ -—{dr,x + Bs~Br)-- —{r,x + Bs-Br)dr ids 
Jo oxi 2 Jo oxi J J 



—^{x + Bs)e'^'ds 
ox: 



ft Qj^ PS Qp Y /"^ 

+ 2 1 —{x + Bs)e^^[ —{dr,x + Bs-Br)-- —{r,x + B,,-Br)dr]ds 



9xi 



dxi 



2 Jq dxi 



Using that 



dxi 

we obtain easily 



f dF ^ , ^ ^ , 1 r dd , „ „ , , 

/ —{dr,x + Bs-Br)-- —(r,x + Bs-Br)dr, 
Jo dxi 2 Jq dxi 



V{t,x) = h{x) + / V{s,x)F{ds,x) + - Y] 

Jo 2 .^^ Jo 



— j{s, x)ds . 



dx. 



This shows that under some the additional regularity conditions on F and h the process u 
defined by (13.51) is a strong solution to Equation (13. 3p . and also a mild solution. 

Step 3. Consider now the case of a general semimartingale F. For any e > we define 



M'{t,x) 



M{t,y)p,{x - y)dy, 



B\t,x) = / B{t,y)p,{x-y)dy, 



and h'^^x) = J^^ h{y)pe{x — y)dy. It is easy to check that is infinitely differentiable 
and it has polynomial growth together with all its partial derivatives. Also F'^{t,x) = 
M'^{t,x) + B'^{t,x) is a C'^-semimartingale with local characteristic given by 



¥{t, x)= hit, y)pe{x - y)dy, 

and ^ 

a^(t,x,y) = / / a(s,x - zi,y - Z2)pe{zi)p^{z2)dzidz2ds. 

Jo JR2d 

It easy to check that a*" and ¥ satisfy the estimates (13. ip and (13. 2p respectively, the partial 
derivatives D^Dya^(t,x,y) satisfy the estimate (13. ip for all multi-index a with 1 < |a| < 2, 
and the functions D'^b'^{t, x) satisfy the estimate (13. 2p for all multi-index a with 1 < |a| < 2. 
From Step 2 it follows that 

V'{t,x) = E^{h'{x + Bt)exp{ [ F'{dr,x + Bt - B,.) -\ f d'{dr,x + Bt - Br))} 

Jo 2 Jq 

is the strong solution to 

dV 1 OF' 

-^{t,x) = -AV%x) + V^^{t,x) ^3_^2) 
V'{x,0) = h'{x) 

As a consequence, it is also a mild solution to (I3.12p . namely, 

V'{t,x)= [ pt{x-z)h'{z)dz+ [ ( [ pt-r{x - z)V'{r,z)F'{dr,z)^ dz. 

jRd jRd \Jo J 

Finally, we are going to take the limit as e tends to zero in each term of the above expression 
in order to deduce the Feynman-Kac formula of V. The estimate (13.10 implies 

sup -Eexp I p / |a'^(r, x + Bt — Br)\dr ] < oo, . 

e>0 \ Jo J 

for all p > 1 and, as a consequence, V^lt, x) converges to V{t, x) in for all p > 1. Clearly, 

lim / pt{x — z)h''{z)dz = I pt{x — z)h{z)dz, 



and 



lim / / pt-r{x — z)V''{r, z)¥{r, z)dzdr = / / pt^ri^ — z)V{r, z)b{r, z)dzdr, 



^-l-O Jo Jr'' Jo Jr'' 

also in for all p > 1. The following limits in are also easy to check: 

lim/" ( [ pt_r{x - z)V'{r,z)[M'{dr,z) - M{dr,z)]) dz = 

and ^ 

lim / (/ pt_r{x - z)[V'{r,z) -V{r,z)]M{dr,z)] dz = 0. 
JRd \Jo J 

This completes the proof of the theorem. ■ 
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4 Stochastic heat equation: Holder continuity of the 
solution 



Consider the following nonlinear stochastic partial differential equation: 
du 1 

— = -Au + b{u) + a{u)W{t,x), t > , x G R"' ^^^^ 
u{0, x) = Uo{x) . 

where W is the Gaussian family introduced in Section 2.1 with covariance function given by 
(12.11) . Let us recall that an adapted random field {u(t,x) ,t > ,x E M.'^} is called a mild 
solution to Equation (14. ip if u satisfies the following integral equation. 

u(t,x)= / pt{x — z)uo{z)dz + / / pt-r{x — z)b{u{r, z))dzdr 
Jr'' Jo Jr'' 

/ pt^r{x - z)(j{u{r,z))W{dr,dz) , (4.2) 
'o Jr<* 

where the stochastic integral is defined as the integral of an H-valued predictable process. 
We are going to impose the following condition on the covariance function. 

(Hi) For each t > 0, 

sup / / pt-six - zi)pt-six - Z2)\qizi, Z2)\dzidz2ds < oo . 

xGiRd Jo JM2<i 

Theorem 4.1 Suppose that b and a are globally Lipschitz continuous functions and suppose 
that the covariance function q satisfies (HI). Let Uq{x) be a bounded function in Mf^. Then 
there exists a unique adapted process u = {u{t,x),t G [0,T],x G W^} satisfying (j-^.^l ). 
Moreover, 

sup E\u{t, x)\P < oo, \/p>2. (4.3) 

te[o,T],xeK'* 

Proof. Fix p > 2. Let Bp be the Banach space of all adapted random fields u such that 
\\u\\p < oo, where = sup E\u{t,x)Y'. On Bp, define the following mapping 

tG[0,T],xeK'* 

'^{u){t,x) := I pt{x — z)uo{z)dz + / / pt-r{x — z)b{u{r, z))dzdr 
Jr^ Jo Jr'^ 

+ / / pt-r{x — z)a{u{r, z))W{dr, dz) . 
Jo Jr-^ 

It is straightforward to obtain 



E\^{u)-'^{v)f{t,x) <C 



E[ I / pt-s{x — z)\u{s, z) — v{s, z)\dzds 

-'R'* 



+ E{ j j pt-s{x - Zi)\u{s,Zi) -v{s,Zi)\pt-s{x - Z2)\u{s,Z2) -V{S,Z2)\ 

'0 Jr^'' 
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X \q{zi, Z2)\dzidz2ds 
Taking the supremum with respect to t and x, we have 

II^(m) - <C f \\u- vf^ds < CT\\u - v\\l . 

Jo 

Consequently, \1/ is a contraction mapping on Bp when T sufficiently small. This proves the 
existence and uniqueness of the solution for some small T. From the above argument it 
is clear that the T such that \E' is a contraction is independent of the initial value of the 
solution. This can be used to to show the existence and uniqueness of the solution for any 
T. The inequality fl4.3p follows in a similar way. ■ 

Now we apply the factorization method to obtain the Holder continuity of u. Fix an 
arbitrary a G (0, 1) and denote 



Y^{r,z)= / pr.siz-y)aiuis,y)){r-s)-''Wids,dy). (4.4) 

Jo JR'i 

The semigroup property of the heat kernel and the stochastic Fubini's theorem yield 

Pt-six - y)a{u{s, y))W{ds, dy) 

Pt-rix ~ z){t - r)°''^Ya{r, z)dzdr. (4.5) 



Jr'^ 
smlvra) '* 







vr 

Consider the following stronger condition on the covariance function. 
(Hla) There exists 7 > — 1 such that for each t > 0, 



sup / pt{x - Zi)pt{x - Z2)\q{zi, Z2)\dzidz2 < Ct^ . 



Lemma 4.2 Let the assumptions of Theorem\JJjbe satisfied. Assume the covariance func- 

1 + 7 

tion q satisfies (Hla). Then for any fixed T > 0,p > 1, a G (0, — - — ), we have 

sup E{\Ya{r, z)\P) < 00. 

r€[0,T],zi 



Proof. Since sup E{\u{r, z)\^) < 00 from Theorem 14.11 and a is Lipschitz continuous, 

rG[0,T],zGRd 

we have 

sup E{\a{u{r, z))\^) < 00. 

re[0,T],zeR'' 

Then we can write 

E\Y^ir,z)f <c(e [ [ p,_,{z - y,)pr^s{z - y2) 

\ Jo JR2d 

X a{u{s,yi))a{u{s,y2))ir - s)~^°'q{yi,y2)dyidy2ds 
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<C sup E{\a{u{r,z))\^) 

re[0,T],zeM<* 

x( / / Vr-s{z -yi)pr~-s{,z -y2){r - s)~'^"\q{yi,y2)\dyidy2ds 
<C (^j^ (r - sf-^'^d^ ' < oo. 



Equation (14. 5p and Lemma 14.21 constitute the main ingredients to prove the following 
theorem concerning the Holder continuity of the solution u. 

Theorem 4.3 Suppose that b and a are globally Lipschitz continuous. Assume (Hla) and 
suppose that uq{x) is bounded and p-Holder continuous. Then the solution u to the equation 
( [^. j[ ) is a.s. Pi-Holder continuous in the time variable t and fi^-Holder continuous in the 

space variable x for any f3i G (0, -[p A (1 + 7)]) and (32 G (0,p A (1 + 7)), respectively. 

Proof. It suffices to follow the idea of the proof of Theorem 2.1 in [15], and we provide a 
sketch of the proof for the reader's convenience. The proof contains two parts for time and 
space variables, respectively. 
Part I 

Fix T,h > and p G [2, 00). First we show that 



sup snp E{\u{t + h,x)-u{t,x)\P) <C{p,T)h'^P, 

0<t<TxGR'' 



(4.6) 

for any i] G (0, |[p A (1 + 7)]). Let be as defined in fl4.4p with a G (0, ^-^) and denote 
Ptf{x) = Jf,dPtix - z)f{z)dz. We have 

4 

E{u{t + h,x) -u{t,x)f) < C{p,a)^Ii{t,h,x), 



1=1 



where 

Ii{t,h,x) = \Pt+hUoix) - PtUo{x)f, 



hit, h,x) = E 
hit, h,x) = E 

hit, h,x) = E 



JM.'i 
t+h 



t+h 



[pt+h-rix -z)it + h- r)"-i - pt.rix -z)it- r)"-i]r,(r, z) 



dzdr J , 



Pt+h-rix - z)it + h- r)" Y^ir, z) 



dzdr 



Pt+h-rix - z)biuir, z))dzdr - I I pt-rix - z)6(u(r, z)) 

JR'* Jo Jw 



dzdr 



For the term hit,h,x), using the fact that mq is p- Holder continuous we have hit,h,x) < 
Ch^ . For any a G (0, i^^) set ■ip"it,x) = By Holder's inequality and Lemma 

we have 



hit,h,x)<c(^j^ 



t r \ P 

\ip"it + h — r,x — z) — ip^it — r,x — z)\dzdr 
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Set 

l2,iit,h,x 

and 

l2,2{t,h,x) = 




exp 



2{t-r) 




exp 



X — z\ 



2{t + h-r) 



exp 



X — z\ 



2{t-r) 



dzdr. 



Then l2(t,h,x) < C(/2,i(t, /i, x)^ + l2^2it,h,xy), and using the same arguments as in the 
proof of Theorem 2.1 in |15j to estimate the terms /2,i(t, h,xy and /2,3(t, h,xy, we obtain 
for G (0, a) that l2{t, h, x) < Ch'^^. By Holder's inequahty and Lemma [4.2^ we have 

hit,h,x) <C (^J^ J pt+h-rix- z){t + h-ry'^dzdr^ <Ch"P. 

A change of variable yields 

h{t,h,x) < C{h^i{t,h,x) + J4,2(t, /i,x)), 



with 



/4 h,x) = E 



h^2{t, h,x) = E 



eh 




Pt+h-r{x — z)b{u{r, z))dzdr 



'0 JR-i 




JR'i 



Pt-r{.x — z) \b{u{r + h,z)) — b{u{r, z))]dzdr 



Since h is Lipschitz, using Holder's inequality and Equation (14.31) . we have 

h,i{t,h,x) < ChF. 
The Lipschitz property of h also implies 



h,2{t,h,x) < / snp E{\u{r + h, z) — u{r, z)\^)dr. 

Jo zi 



Putting together all estimations for Jj, z = 1, . . . , 4, we obtain 

sup E{\u{t + h,x) -u{t,x)\P) < C ('/iP™"(f + /" sup E{\u{r^h,x)-u{r,x)\^)dr 

Since < < a < i^^, the estimate (14. 6p follows by Gronwall's Lemma. 
Part II 

Now consider the increments in the space variable. We want to show that for any T > 0,p G 

[2, cx)), X, a G M'^ and G (0, p A (1 + 7)), 



sup sup E(|M(t,x + a) -M(t,x)|P) < Cjal^'P. 

0<t<T xi 



(4.7) 
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Fix a e (0, Ee have 



E{\u(t, X + a) — u(t, x)\^) < C Ji{t, X, a), 



with 



Ji{t,x,a) = \PtUo{x + a) - PtUo{x)\^ 



J2(t, x,a) = E 
J3{t,x,a) = E 



Jw^ 

t 



[ip'^{t — r, X + a — z) — ^"{t — r, X — z)]Ya{r, z)dzdr 



~'K'' 



[pt-r{x + a — z) — pt-r{x — z)]b{u{r, y))dzdr 



It is easy to show that Ji(t, x, a) < C\a\f^. For the term J2{ti x, a), first we have using the 
mean value theorem, 



\i;''{t-r,x + a- z)-ip''{t-r,x- z)\dz < C{t-r) 



where rj G (0,1). Again by Holder's inequality and Lemma W?2\ for a G (0, -|^),r7 G 
(0, 2a A 1), we deduce 



J2(t, X, a) < C 



\i)''{t - r,x + a - z) - ij"{t - r,x - z)\dzdr^ < C\a\^'P . 



Finally, by a change of variable, the Lipschitz property of b, and Holder's inequality, 



J^{t,x,a) < E 



Pt-r{x — y) [6(M(r, z + a) — b{u{r, z)]dzdr 



< C snp E{\u{r, a + z) — u{r, z)\^)dr. 



2 



Then fl4.7p follows from the Gronwall's lemma and the estimates of Jj, i = 1,2, 3. 
Here are two examples. 



Example 4.4 A similar Holder continuity result was obtained in [T5\l in the case of an 
homogeneous covariance function q{x,y) = q{x—y) > 0, where q is a nonnegative continuous 
function on M'^\{0} such that it is the Fourier transform of a non-negative definite tempered 
measure fi on M"^, and for some rj G (0, 1) we have 



< oo. 



(l + I^P)^ 

This condition implies (Hla). In fact, we can write 

Pt{x - zi)pt{x - Z2)qizi - Z2)dzidz2 = / pt{x - z - Z2)pt{x - Z2)q{z)dz2dz 



P2t{z)q{z)dz 



< 



-2t(52 + l) 



^{di) < Ct~^ 



+ 



2^1 r, • 
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Theorem 14.31 can be applied to noises which do not have an homogeneous spatial covari- 
ance like the following example. 

Example 4.5 Consider the case where d = 1 and the covariance structure in space is that 
of a bifractional Brownian motion with parameters if G (0, 1), G (0, 1], that is, 

where 2KH > 1. Then, B^'^{t,x) = M^(l[o,t]x[o,x])j is a bifractional Brownian motion in 
X > for each fixed t, and formally, W{t, x) = -^B^'^(t, x). Then 

\qix, y)\ < C[\xr^-' + \yr^-' + \x - y^^-'] 

and 7 = HK — 1 G (—1, 0). Thus, Theorem \4.3\ can be applied to this case. 

5 Stochastic heat equation: Regularity of the density 
of the solution 

In this section we consider again the solution u{t, x) to (14. ip . and we will impose the following 
condition on the covariance q{x,y). 

(H2) q is 7o-Holder continuous for some 70 > 0, and for some (3 G [0, 2) 

\q{xi,X2)\<C{l + \x,\^ + \x2f). 

In this case we can assume that the random field {VFj has a density with respect to the 
Lebesgue measure on M.'^. That means, we suppose that there exists a zero mean Gaussian 
random field {Wi{t,x),t >0,xE W^} with covariance 

E{Wi{t,x)Wi{s,y)) = {sAt)q{x,y), 

such that Wt{ip) = J^d ^{x)Wi(t, x)dx, for any ip G Co(ffi'^), where q{x, y) is positive definite, 
namely, f^2d(l{x,y)f{x)f{y)dxdy > for all / G L^(M.'^,dx). The additional regularity 
conditions imposed on q have allowed us to introduce the density process Wi(t,x), which is 
a Brownian motion in the time variable and it has the spacial covariance q. 

From a Theorem of Mercer's type (section 98 on page 245 in [H]) we know that if 
J^2d \q{x,y)\'^dxdy < 00, then 



n=l 



where e„, n = 1, 2, ■ ■ ■ is an orthonormal basis of L^(]R'^, dx) and Yl'^=i < The positive 
definite property of q{x,y) implies A„ > 0. If we take C{x,y) = V^^n(^)^n(l/), 
then q{x,y) = f^dc{$,,x)c{^,y)d^. Thus it is without loss of generality for us to assume 
that q{yi,y2) = J-^dC{^,yi)c{.^,y2)d^ for some c{^,y). Furthermore, we assume c{x,y) has 
polynomial growth. 

The following is the main result of this section. 
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Theorem 5.1 Assume that q is 'jQ-Holder continuous for some 70 > and satisfies (Hla). 
Suppose 

g(xi,X2) <C(l + |xi|'^ + |x2|^) /or some /3 e [0, 2) . (5.1) 

Let uq be bounded and p-Holder continuous for some p > 0. Suppose that there is a xq E 
such that q{xo,Xo) > and o"(uo(xo)) 7^ 0. Then, 

(1) Ifb and a are continuous differentiahle functions with bounded first order derivatives, 
for any t > and x G W^, the probability law of u{t,x) is absolutely continuous with 
respect to the Lebesgue measure. 

(2) Ifb and a be infinitely difjerentiable with bounded derivatives of all orders, then for 
any t > and x G M*^, the probability law ofu{t,x) has a smooth density with respect 
to Lebesgue measure 

Proof. First we claim that for aU {t,x) the random variable u{t,x) belongs to the 
Sobolev space D^'^ under condition (1), and to the space D°° under condition (2). This 
follows from standard arguments and we omit the proof (see, for instance ^IJJ, Proposition 
2.4.4 in the case of the stochastic heat equation). On the other hand, the Malliavin derivative 
Ds,yu{t, x) satisfies the linear stochastic evolution equation 

Ds,yu{t,x) =pt-s{x - y)a{u{s,y)) + / / pt-r{x - z)b'{u{r, z))Ds,yu{r, z)drdz 

jR-i Jo 

+ / pt-r{x — z)a'{u{r, z))Ds^yu{r, z)Wi{dr, z)dz . 
Jr'' Jo 

Denote 

F := \\Du{t,x)\\l^^ = / D,^y^u{t,x)Ds,y.^u{t,x)q{yi,y2)dyidy2ds , 

Jo JM2d 

where || ■ Wy,^ is the Hilbert norm introduced in Section 2. We are going to show only the 
statement (2), and the first one follows from similar arguments. It suffices to show that 
E [F~^] < 00 for any p > 1. We divide the proof into two steps. 

Step 1. Introduce Vs,({t,x) = / c{^,y)Ds^yu{t, x)dy. Then we can write 

Jw 

F = [ [ Vs^^{t,xfdids. 
Jo Jw 

For any fixed (s,^), the random field {Vs^^{t,x),t > s, x G M°'} satisfies the following linear 
stochastic heat equation for t > s, and x G M°', 

= l^Vs,, + b'{u)Vs,, + a'{u)V,,,^{t, x), 
ys,ds,x) = c{^,x)a{u{s,x)) . 

Consider the continuous semimartingale {F{t,x),t > 0, x G M'^} given by 

F{t,x) = / b'{u{r, x))dr + / a'{u{r, x))Wi{dr, x). 
Jo Jo 
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The local characteristic of this semimartingale are b{t,x) = b'{u{t,x)) and 



a{t,x,y)= / a'{u{r,x))a'{u{r,y))q{x,y)dr . 



Notice that conditions (13. ip and (13. 2p hold because b' and a' are bounded and q satisfies 
(15. ip . Then, Theorem 3.1 gives an explicit Feynman-Kac formula for the above equation. 
This means that we have 



V,^^{t,x) =E 



B 



c(^, X + Bt-s)(y{u{s, X + Bt-s)) 



X exp { I F{dr, x + Bt^s - Brs) ^\ a(r + x, Bts - Br-s)dr 



Step 2. Let 

t;B) = F{dr, x + Bts - Brs) - ^ d{r + x, Bts - Brs)dr. 



Then 



\Vs,^it,x)\'d^ds= / / E 
Jm.'^ Jo Jm.'' 



B,B 



c{^,x + Bts)c{^,x + B 



X a{u{s, X + Bt-s))o{u{s, x + Bts)) exp{F(s, t; B) + F(s, t] B)} 



dC,ds 



E 



B,B 



q{x + Bts,x + Bts 



ds 



X a{u{s, X + Bts))criu{s, x + Bts)) exp{y(s, t; B) + Y{s, t; B)} 



H{s)ds, 



where is a standard Brownian motion independent of -B. If we can show that E{H{0) ^) < 
oo for all p > 1, and H{s) is Holder continuous, then by Lemma [5^ below we deduce 



ft p \-p 

E{ I / \Vs,^{t,x)\^d^ds] =E 
'o 



H{s)ds < oo 



for all p > 1. The Holder continuity of H{s) can be verified from the following inequality: 



E\H{s^)-His2)f<C\s2-s, 



I £ mm{p,7o,l+7} 



where C is determined by 



sup \E\q{x + Bts, X + Bts)f'', E\a{u{s, x + Bts))f'' , E exp{8pY{s, t-B)}]. 
se[o,t] I J 

It remains to show that E{H{0)~p) < oo. Notice that 

H{0) = E^-^ (g^ exp{y(0, t; B) + Y{0, t; B)} 



where 

Gx = q{x + Bt, X + Bt)a{uo{x + Bt))a{uo{x + Bt)). 
We can write, by Jensen's inequality, 



--E 



E 



B,B 









< 




E 



G^exp{Y{0,t;B) + Y{0,t;B)} 
\G,\ sign(G',) exp{r (0, t; B) + F (0, t; B)] 

G,\ exp{-p(r (0, t; B) + F (0, t; E))} 



Our nondegeneracy hypotheses imply that E^'^G > 0, and this allows us to conclude the 
proof. ■ 

Lemma 5.2 Let {St,0 < t < 1} be a non-negative stochastic process. If ESq"' < oo for 
some a > 0, and sup — 5'o| < Gf^ where G is a positive random variable with EG^ < oo 

0<s<t 

for some b > 0, then we have 



E 



Stdt 



< oo, for < p < ab^l (a + 6 + 67) . 
In particular, if a and b can be arbitrarily large, then p can also be chosen arbitrarily large. 
Proof. Let a, ^ > 0, where a + P < 1 and 6^7 — ba > aa, and < e < 2"+^^"^. We have 









\ r' 1 




p 


[ Stdt < e 


< P 


/ Stdt < e, -So > 


+ P[Sq< e'^ 




Jo 




Jo 





<p 



sup \St - Sq\ > ie" 



P [So" > e"""] 



<2^€-^''E I sup l^t - 5'o|'' I + e'^'ESo" < C (e^'^^-^'' + e"") < Ce" 

, 0<t<el^ 



Then E 
such that 



lo Stdt 



< 00, for < p < aa. The lemma follows with the choice of a and /3 



a < 67/(0 + 6 + 67) and /3 = (a + 6)/(a + 6 + 67). ■ 
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